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ABSTRACT: In a majority of mammals, male infants have heavier body mass and grow faster than female infants. Accordingly, male 
offspring nursing requires a much greater maternal energy contribution to lactation. It is possible that the maternal-fetal 
immunoendocrine dialog plays an important role in female preparation for lactation during pregnancy. Immune system genes are an 
integral part of gene regulatory networks in lactation and tumor necrosis factor alpha (TNFa) is a proinflammatory cytokine that also 
plays an important role in normal mammary gland development. The aim of this study was to evaluate the influence of the sex of calf 
and/or the -824A/G polymorphism in the promoter region of TNFa gene on milk performance traits in Black Pied cattle over the course 
of lactation. We also studied the allele frequency differences of -824A/G variants across several cattle breeds, which were bred in 
different climatic conditions. The G allele frequency decreased gradually over the course of lactation events in the Black Pied dairy 
cattle because of a higher culling rate of cows with the G/G genotype (p<0.001). In contrast to the genotypes A/A and A/Q cows with 
G/G genotype showed significant variability of milk and milk fat yield subject to sex of delivered calf. Milk yield and milk fat yield 
were significandy higher in the case of birth of a bull calf than with a heifer calf (p<0.03). The G allele frequency varies from 48% to 
58% in Grey Ukrainian and Black Pied cattle to 77% in aboriginal Yakut cattle. Our results suggest that the TNFa -824A/G gene 
polymorphism may have an influence on the reproductive efforts of cows over the course of lactation events depending on the sex of 
progeny. Allocation of resources according to sex of the calf allows optimizing the energy cost of lactation. This may be a probable 
reason for high G allele frequency in Yakut cattle breeding in extreme environmental conditions. Similarly, the dramatic fall in milk 
production after birth of a heifer calf increases the probability of culling for the cows with the G/G genotype in animal husbandry. (Key 
Words: Cattle, TNFa Gene, Single Nucleotide Polymorphism, Milk Production, Sex of Progeny) 



INTRODUCTION 

Sexual size dimorphism becomes apparent in heavier 
body mass of males in most mammal species (Lark et al., 
2006). These differences are formed during prenatal 
development and increase during the ontogenesis (Loos et 
al., 2001). Accordingly, breastfeeding of male offspring 
requires a higher maternal investment in lactation. 
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Physiological lactation support depending on the sex of 
newborns occurs during pregnancy (Galbarczyk, 2011). 
Later the differences in milk productivity are caused by 
stronger breast stimulation of male offspring (Glutton- 
Brock et al., 1981). Maternal organism receives information 
about the fetus male sex through the maternal iimnune 
response to HY fetal antigens (Nielsen, 2011). Intrauterine 
development is also regulated by different cytokines. Tumor 
necrosis factor alpha (TNFa) plays an important role in 
development of mammary gland (Warren et al., 2009), and 
inflammatory response of uterus after mating (Sanford et al., 
1992). TNFa is abortifacient in high doses (Chaouat et al., 
2007), but is necessary in low doses for the local tissue 
remodelling required for implantation (Chaouat et al., 2004). 
The possible function of TNFa may be to boost death 
signaling in order to kill the embryo if initial damages 
triggered by detrimental stimuli may culminate in structural 
anomalies, and stimulate protective mechanisms if the 
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repair of these damages may prevent maldevelopment 
(Toder et al., 2003). The variabiUty of milk traits according 
to sex of offspring can facilitate adaptive evolution in native 
mammals, but this variability may be a harmful factor for 
dairy cattle breeding. The aim of this study was to evaluate 
the influence of the sex of calf and the -824A/G TNFa 
genotypes on milk traits in Black Pied cattle and the allele 
frequency differences across domestic cattle populations. It 
was hypothesized that sex or genotype could influence milk 
production and selection can shift the allele frequencies. 

MATERIALS AND METHODS 

Animals 

Blood samples of Black Pied breed (n = 107) were 
collected in the experimental farm "Elitnoe" of the Siberian 
Research Institute of Animal Husbandry SD RAAS 
(Novosibirsk region). The data on milk performance traits 
for each Black Pied cow were taken from breeding 
documentation. On the experimental farm "Elitnoe", 
lactating dairy cows were housed either in a free-stall barn 
or a pasture. The cows were milked twice a day. Automatic 
milking machines allowed obtaining milk yield data for 
each cow. The milk fat content was determined once a 
month. The milk yield, milk fat yield, and fat percentage 
were calculated by the A4/2 test interval method for the 
standard 305-d lactation length (ICAR, 2012). The calves 
were separated from their mothers within a day after their 
birth. After separation, young dairy calves subsisted on 
liquid milk, a commercial milk-replacer, or colostrum. 

Blood samples of Grey Ukrainian breed (n = 48) and 
Yakut breed (n = 33) were collected from an experimental 
agricultural farm of SB RAS (Altai Repubhc). Grey 
Ukrainian and Yakut cattle herds were raised at this 
experimental farm according to the conservation program of 
livestock genetic resources and were not used in nulk 
production. 

Genotyping of DNA variants 

DNA was isolated by the standard methods using 
proteolytic treatment followed by phenol extraction. A 
fragment of 5 '-flanking region of the bovine TNFa gene 
between positions 3090 and 3234 (GenBank accession no. 
Z14137) was amplified from genomic DNA by PGR with 
forward primer 5'-CCGAGAAATGGGACAACCT-3' and 
reverse primer 5'-GCCATGTATCCCCAAAGAAT-3'. The 
amplicon was digested with EcolCR 1 (SibEnzyme, Russia) 
and subsequently submitted to poly aery 1 amide gel analysis 
for the restriction fragment length polymorphism 
determination of the -824A/G site. Genotype nomenclature 
refers to the respective SNP in position 3169 (GenBank 
accession no. Z14137). Fragment lengths were 81/64 bp for 
A/A homozygotes, 145 bp for G/G homozygotes, and 



145/81/64 bp for A/G heterozygotes. 

Statistical analysis 

Allele frequencies were calculated by direct counting. 
Allele frequencies were compared using contingency table 
(2x2) chi-square tests. Average numbers of lactations 
between genotypes were compared using nonparametric 
Kruskal-WaUis test. Associations of the SNP polymorphism 
-824A/G in the promoter region of the bovine TNFa gene 
with milk performance traits were evaluated with analysis 
of covariance (ANCOVA, covarying for number of lactation 
episodes) using genotype and sex of the calf as independent 
variables. Post hoc analysis was performed with correction 
for multiple comparisons by Fisher's LSD test. Taking into 
consideration that after the fourth lactation there was less 
than 20% of initial livestock remaining, the milk 
performance traits were analyzed using the data from first 
to fourth lactations only. Since one cow could have several 
lactations in its lifetime, the number of genotypes exceeded 
the number of animals in the statistical analysis. The data 
were analyzed with the GLM procedure of STATISTICA 
8.0 software. 

RESULTS 

No significant differences were observed in frequency 
of G allele of the TNFa -824A/G gene polymorphism 
between Grey Ukrainian (47.9%) and Black Pied (57.9%) 
cattle breeds. At the same time G allele frequency was 
increased to 77.3% in Yakut cattle breed. The G allele 
frequency decreased gradually from 57.9% in primiparous 
heifers to 38.0% in fifth lactation and older cows of the 
Black Pied dairy herd which is likely due to differential 
culling of herd at the farm "Elitnoe" (Figure 1). This 





0,7 1 






0,6 ■ 




0,5 ■ 




6 


0,4 ■ 






o 


0,3 ■ 


t' 


c 


0,2 ■ 


i) 


3 

cr 


0,1 ■ 


u 


u 


0,0 ■ 



1 2 3 4 >5 
Number of lactation 

Figure 1. Temporal changes of the G allele frequency of the 

TNFa -824A/G polymorphism per number of lactations in Black 
Pied cattle which is likely due to differential culling of herd, n = 
31-107 animals per point. Significant differences in the G allele 
frequencies were observed between the first and fourth lactations 
(x^ = 4.02, p<0.045) and between the first lactation and combining 
data for fifth and following lactations {x'= 10.85, p<0.02) as well. 
The G allele frequency was negatively correlated with number of 
lactation (rs = -0.99, p<0.002). 
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decrease was due to frequency of the G/G genotype 
decreasing faster than A/A and A/G genotypes (Figure 2). 
The average number of lactations during Hfetime in cows 
with the A/A, A/G and G/G genotypes were 2.92±0.21, 
2.96±0.18, and 1.99±0.13, respectively. There was a highly 
significant difference among the three genotypes in this 
parameter (Kruskal-Wallis test H2,3i3= 14.89, p<0.001). The 
percentages of bull calves to heifer calves in cows with the 
A/A, A/G and G/G genotypes were 47.1%, 43.1%, and 
45.3%, respectively. 

The 305-d milk yield and milk fat yield increased from 
first to fourth lactation in cows of all genotypes (Figure 3A 
to C). After fourth lactation, either a tendency to decrease or 
no further changes in milk yield and milk fat yield were 
observed. The percentage of fat didn't change significantly 
throughout the entire life spans of the cows. ANCOVA with 
niunber of lactation as a covariate was used to evaluate the 
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Figure 2. Temporal changes in the percent of total remaining for 
each TNFa -824A/G genotype per number of lactations in Black 
Pied cattle which is likely due to differential culling of herd, n = 2- 
46 animals per point. The standard error of the mean for genotype 
trequencies is displayed. The frequency for each genotype in 
primiparous heifers is taken as 100 percent. 
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Figure 3. Data on milk performance traits in cows with different TNFa -824A/G genotypes. Data are reported as mean+SEM. A, B, C: 
temporal changes of 305-d milk yield, 305-d percentage of fat, and 305-d milk fat yield, respectively; D, E, F: ANCOVA results for the 
same traits. Number of lactations is used as a covariate and data from first to fourth lactation are combined for ANCOVA. * p<0.02 
relative to the G/G group corrected for multiple comparisons using Fisher's post-hoc LSD test. 
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Table 1. Influence of genotype and sex of calf on milk performance traits in Black Pied cattle breed according to ANCOVA with number 
of lactation as a covariate 



Factor 


305-d milk yield 


305-d percentage of fat 


305-d milk fat yield 


Genotype 


F2,255 = 4.63,p<0.01 


^2,255 = 0.03, p<0.97 


^'2,255 =4.45,p<0.01 


Sex of calf 


Fi,255 = 2.51,p<0.11 


^'i,255 = 0.08, p<0.77 


Fi,255 = 2.27,p<0.13 


Intetraction genotypexsex 


^'2,255 = 5.07, p<0.007 


^'2,255 = 1.99, p<0.14 


F2,255 = 4.70,p<0.01 



influence of sex of the calf and TNFa genotype on milk 
performance traits. The influence of a cow's genotype and 
genotypexsex of calf interaction on milk yield and milk fat 
yield has been established (Table 1). Milk yield and milk fat 
yield were significantly lower in heterozygotes for the 
TNFa -824A/G gene polymorphism compared to G/G 
homozygotes (p<0.02) and had a tendency to decrease 
compared to A/A homozygotes (p<0.06) (Figure 3D to G). 
At the same time sex of the calf alone didn't influence milk 



performance traits. But it had significant effect on these 
traits in interaction with cow's genotype (Table 1). It turned 
out that the milk yield and milk fat yield in cows with the 
G/G genotype were significantly higher in case of birth of a 
bull calf than that in case of birth of a heifer calf after 
correction for multiple comparisons (Figure 4A to C). Sex 
of the calf did not have an influence on milk performance 
traits in cows with the genotypes A/A and A/G. 

Sex of the calf could influence the milk performance 




Figure 4. Influence of sex of the calf on the milk performance traits of cows with different TNFa -824A/G genotypes. Average data from 
first to fourth lactation are shown. Data are reported as mean±SEM. A, B, C: ANCOVA results for the 305-d milk yield, 305-d percentage 
of fat, and 305-d milk fat yield, respectively. Number of lactations is used as a covariate. D, E, F: ANCOVA analysis of residual 
variances of the same traits. * p<0.00004, ** p<0.02 relative to the buU calves corrected for multiple comparisons using Fisher's post-hoc 
LSD test. 
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traits because cows that usually give birth to bull calves had 
higher milk yield and milk fat yield than cows that usually 
give birth to heifer calves. To answer the question of 
whether the cow's milk performance traits increase or 
decrease after giving birth to a bull or heifer calves, analysis 
of residual variances of milk performance data was 
performed with centering the variables around its means for 
the entire life spans of individual cows. Only data from the 
cows that had two and more lactations and delivered both 
bull and heifer calves were centered. 

ANCOVA with number of lactation as a covariate 
showed that residual milk yields were significantly 
influenced by sex of the calf and genotype x sex interaction 
(^"1,133 = 4.78, p<0.03 and ^2,133 = 2.11, p<0.05, respectively). 
According to post hoc analysis, cows with the G/G 
genotype showed significantly higher milk yield in case of 
birth of a bull calf than that in case of birth of a heifer calf 
after correction for multiple comparisons (+447.4±152.8 kg 
and -362.2±205.5 kg, respectively) (Figure 4D). Sex of the 
calf and TNFa genotype didn't influence the percentage of 
fat. Residual milk fat yield, as weU as residual milk yield, 
were significantly influenced by sex of the calf and 
genotype x sex interaction (^1133 = 5.29, p<0.02 and F2J33 
= 3.47, p<0.03, respectively). Cows with the G/G genotype 
also showed significantly higher milk fat yield in case of 
birth of a bull calf than that in case of birth of a heifer calf 
after correction for multiple comparisons (+19.3±6.2 kg and 
-15.6±7.8 kg, respectively) (Figure 4F). 

DISCUSSION 

The TNFa gene is present as a single copy gene 2,773 
base pairs (bp) in length on bovine chromosome 23. The 
gene consists of four exons and three introns. mRNA 
product length is 1,713 bp. There are 28 single nucleotide 
polymorphisms (SNP) of bovine TNFa gene in the dbSNP 
database, but validation status for most of them is unknown 
(Sherry et al., 2001). Two SNP polymorphisms of 
TNFalpha gene in the exon and promoter regions were 
found to have a strong association with the early first 
ovulation after parturition in the high-producing dairy cows 
(Shirasuna et al., 2011). Since no effect of these 
polymorphisms on the gene function has not yet been 
proved, it may be possible that these markers are not 
functional variants, but may simply be in linkage 
disequilibrium with the true casual polymorphism(s). So far 
only one SNP (-824A/G) of TNFa gene was found may 
affect the functionality of the gene, since it influences the 
transcriptional promoter activity in vitro (Konnai et al., 
2006). Moreover, the frequency of -824G allele was higher 
in individuals with bovine leukemia virus (BLV) induced 
lymphoma than in asymptomatic carrier individuals. Thus 
-824A/G SNP has been genotyped in our study. 



We studied the frequency of this SNP among several 
cattle breeds, which were bred in different cUmatic 
conditions. The Yakut cattle live above the Arctic Circle at 
all times (Granberg et al., 2009). Yakut cows have low 
milkabiUty but high resistance to the hardest Arctic 
conditions, where winter temperature drops up to -60°C and 
summer temperatures reaches +40°C. During 10 months 
temperatures remain below zero. In winter the cattle are 
kept in unheated cowsheds and even at the temperature 
-50°C, the cattle are sent to drink from ice-holes. Grey 
Ukrainian cattle initially were a breed meant for producing 
meat and working the fields: it could work for 12 h straight 
in the fields. In 1917, it was recognized as the best beef 
breed in Russia. Currently, around 100 head of these cattle 
are conserved in Siberia. Ecophysiological and molecular 
phylogenetic characteristics also indicate a close 
relationship of aboriginal Yakut and Grey Ukrainian cattle 
to the aurochs (Semenova et al., 2000; Kantanen et al., 
2009). Before domestication aurochs were widespread in 
Eurasia from the Scandinavian peninsula in the north to 
Northern Africa in the south (van Vuure 2005), and 
probably they were well adapted to the harsh enviroimient. 
Black Pied breed is the major breed of dairy cattle utiUzed 
in milk production in the Russia. 

Our study of Russian cattle breeds shows that G allele 
frequencies of the TNFa -824A/G gene polymorphism in 
Grey Ukrainian and Black Pied breeds are 47.9% and 
57.9% respectively. These values are in good agreement 
with previously reported value of 51.9% for Japanese black 
and Holstein-Friesian cattle breeds (Konnai et al., 2006). 
Difference in G allele frequency between the breeds can be 
explained by genetic drift and founder effect, because the 
effective population size during breed development usually 
is small. The aboriginal Yakut cattle breed in our study had 
high frequency of G allele - 77.0%. Yakut cattle is well- 
adapted to extreme climates and closely related to the wild 
ancestor of domestic cattle - the aurochs {Bos primigenius) 
(Kantanen et al., 2009). 

The sex ratios of bull calves to heifer calves in cows 
with the A/ A, A/G and G/G genotypes were 47.1%, 43.1%, 
and 45.3%, respectively. These ratios did not differ 
significantly from the overall secondary ratios (at birth) nor 
in Bos taurus (51.52%) (Asdell, 1964), nor in Bison 
bonasus (49.2%) (Jaczewski, 1958). Analyzing age specific 
variation of the TNFa -824A/G gene polymorphism, we 
have found that the G allele frequency decreased gradually 
over the course of lactation events. In addition, the average 
numbers of lifetime lactations in cows with the G/G 
genotype were one-and-a-half times lower than that in cows 
with the genotypes A/A and A/G. This is possible only if 
cows with the G/G genotype have a higher culling rate. But 
what is the cause of such high culling rate, which leads to 
decrease of G allele frequency in dairy herds compared with 
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domestic low yielding cows such as Yakut cattle? 

Although domestic cows can Uve for 20 yrs, those 
raised for dairy rarely live that long, as the average cow is 
removed from the dairy herd around age four and marketed 
for beef (Ng, 2001; Wallace, 2002). On most dairy farms, 
annual culling rates exceed 30% of the herd. These animals 
are culled mainly due to low milk production or infectious 
diseases (mastitis). As to susceptibihty to infectious 
diseases, it was reported that the frequency of the G allele 
of the TNFa -824A/G gene polymorphism was higher in 
individuals with BLV-induced lymphoma than in 
asymptomatic carrier individuals (Konnai et al., 2006). 
Similarly, these authors observed a tendency for increased 
provirus load in cattle with G/G homozygote compared to 
A/A homozygote or A/G heterozygote. These data suggest 
that the TNFa -824A/G gene polymorphism could at least 
in part contribute to the progression of lymphoma in BLV 
infection. BLV is one of the most common infectious 
viruses of cattle with a worldwide distribution. BLV 
infected cows have a higher culling rate due to higher 
susceptibility to opportunistic and infectious agents 
(mastitis, lameness diseases, ringworm infection etc.), 
decreased productivity, and lower reproductivity (Rodriguez 
etal., 2011). 

As to low milk production, it is necessary to note that if 
the cow's yearly estimated milk production level is 20% or 
more below the herd's average it may be enough 
justification to base a culhng decision (Small Farm 
Resource, 2002). According to our data, cows with the G/G 
genotype are characterized by considerable variation of the 
milk yield and milk fat yield between lactations in contrast 
to A/A homozygote or A/G heterozygote, in which sex of 
the calf didn't influence on milk performance traits (Figure 
4). The milk yield and milk fat yield of cows with the G/G 
genotype were decreased by almost 17% in case of birth of 
a heifer calf than that in case of birth of a bull calf. It is 
clear that such a decrease in milk production after birth of a 
heifer calf increases the probabihty of culhng of the cow. 

The reason for variation in milk yield and milk fat yield 
that depends on the sex of calves may be a kind of dialogue 
that goes on during pregnancy between the mother and her 
unborn baby. This dialogue prepares the mother for 
lactation that is adequate for the physiological needs of 
newborn calves. Genotype-restricted influences of sex on 
milk performance traits take place only in cows with the 
G/G genotype. Our analysis of residual variances of milk 
yield, percentage of fat, and milk fat yield over the entire 
life span of individual cows confirmed the contribution of 
fetal sex to variation in the milk production traits. 
According to the analysis of residual variances, lactations 
following the birth of a heifer calves were characterized by 
lower milk productivity in comparison with the 
corresponding individual average value for the entire life 



span of the cow. Similarly, the milk production following 
the birth of a bull calf exceeds the individual's Ufetime 
mean. These differences in lactation in response to the sex 
of newborn calf relate only to milk and milk fat yield. Sex 
of the calf didn't influence the percentage of milk fat. Thus, 
analysis of residual variances of milk performance data 
shows that the sex of progeny significantly influences 
preparation of the mother for lactation. 

In a majority of mammals, male infants grow faster and 
demand increasing milk production from their mothers 
(Lark et al., 2006). It is possible that a male fetus prepares 
the mother's body for more intensive lactation already 
during intrauterine development. As this effect was 
observed in cows with the G/G genotype only, it seems 
most likely that signal transduction in the "mother-fetus" 
system occurs by HY antigens-induced mother's immune 
responses. It is known that several widely expressed genes 
on the Y chromosome code for male-specific proteins that 
provoke an immune reaction in mothers (Graves, 2010). 
TNFa is known to induce a normal mammary gland 
development, but at the same time it is also able to induce 
apoptotic cell death (Locksley et al., 2001). The mean 
transcriptional promoter activity is significantly decreased 
in -824G promoter compared to -824A promoter of the 
cattle TNFa gene in vitro (Konnai et al., 2006). In these 
terms, probable low TNFa production in cows with the G/G 
genotype may enhance the number of mammary epithelial 
cells. Finally decreased TNFa production may also be 
involved in the modulation of Thl/rh2 balance during 
materneil immune response to HY antigens (Herman et al., 
2011). 

Allocation of resources according to needs of male and 
female progeny in cows with the G/G genotype may 
improve reproductive fitness of wild cattle species. This 
idea is supported indirectly by its high frequency in the 
Yakut cattle, Uving in extreme Arctic conditions with 
extremely cold winter and high ectoparasite (mosquito etc.) 
burden in sunmier. 

The efficiency of infant suckUng and milk removal 
governs the milk volume produced in the udder. 
Observations show that male calves suck from their 
mothers more frequently than female calves (Clutton-Brock 
et al., 1981). But in dairy farms calves are separated from 
their mothers immediately after birth and are fed with 
artificial food (More et al., 2010). In particular, this practice 
takes place in the experimental farm "Elitnoe" where our 
study was carried out. Other authors found that the sex of 
offspring had no effect on milk yield (Wohlt et al., 1984). 
But they didn't take into consideration possible influence of 
genetic factors and calf nursing behavior on this trait. 

Another reason for different milk production among 
G/G individuals that is correlated with the sex of the calf 
may be the stable individual differences in milk 



Yudinetal. (2013) Asian Austr alas. J.Anim. Sci. 26:1379-1387 



1385 



performance traits among cows. For example, more 
productive cows may more often give birth to bull calves. 
Actually in humans, the mothers of boys are heavier in 
pregnancy than the mothers of girls and have a higher body 
mass index (Alwasel et al., 2011). However the analysis of 
residual variances of milk performance data with centering 
the variables around its means for the entire hfe spans of 
individual cows disprove this hypothesis. 

Dramatic changes in cell composition and fimction 
occur in the mammary gland during a pregnancy-lactation- 
involution cycle. Immune system genes are an integral part 
of gene regulatory networks in lactation. Milk contains a 
variety of innate defense factors and cellular components 
that provide both protection against infections and allergies 
and it also stimulates the appropriate development of the 
infant's own immune system. Immune genes expressed in 
mammary tissue during pregnancy appear to participate in 
mammary gland development (Faucon et al., 2009; Watson, 
2009). Some genes that have been annotated with respect to 
positive regulation of lymphocyte proliferation and 
cytokine biosynthesis are important to the development of 
the mammary gland during late pregnancy. However, genes 
associated with such terms as 'antigen processing' or 
'antigen presentation' are statistically down-regulated during 
lactation relative to late pregnancy (Lemay et al., 2007). 
Postlactational regression (phagocytosis of cell debris and 
milk droplets) is associated with a phase of acute 
inflammatory response, in addition to cell death (Atabai et 
al., 2007; Watson, 2009). This is consistent with the original 
idea that the mammary gland evolved from the innate 
immune system (Vorbach et al., 2006). 

An association of immune system genes with 
reproductive traits of livestock animals may be considered 
in two respects. First, there may be a positive correlation 
between immune responses and these traits, because of 
animal breeding in the presence of pathogens that would 
favor improving reproductive traits mostly in the 
genetically resistant animals. Secondly, the correlation may 
be negative if a trade-off between immune function and 
reproductive efforts is involved. It is obvious that 
developing and maintaining an effective immune system 
under parasite challenge may be costly and result in reduced 
secondary sexual (reproductive) traits (van der Most et al., 
2011). Unfortunately, there are only few studies analyzing 
the associations of polymorphisms in inamune genes with 
reproductive traits in livestock animals. For example, the 
association of lymphoid enhancer-binding factor- 1 (LEFl) 
gene with the number of functional and inverted teats in 
pigs was described (Tetzlaff et al., 2009). The association of 
the proliferation associated nuclear element 1 (PANEl) gene 
encoding one of the minor histocompatibihty antigens with 
the hematological immune traits influences the birth weight 
in pigs (Huang et al., 2010). The association of SNPs in the 



chemokine CCL2 and caspase recruitment domain CARD15 
genes with milk production traits was found in Canadian 
Holstein cattle (Leyva-Baca et al., 2007; Pant et al., 2007). 
Other authors reported the associations of SNPs in Toll-hke 
Receptor 4 (TLR-4) gene, chemokine receptor 1 (CXCRl) 
gene, CD14 gene, and serum protease inhibitor 
(SERPINAl) gene with milk production in dairy cattle 
(Beecher et al., 2010). This conceptual framework clearly 
warrants further investigation. 

We also observed the low milk yield and milk fat yield 
in A/G heterozygotes compared to both A/A and G/G 
homozygotes, i.e. so-called underdominance effect. Other 
authors described a similar underdominance for some 
reproductive traits in cattle (Lin et al., 1987). The cause of 
this phenomenon is unclear, but it may suppose that reduced 
milk production in heterozygotes results in higher relative 
fitness (heterozygous advantage), possibly increasing 
siu'vival rates in harsh environment. 

In conclusion, we found evidence that the TNFa - 
824A/G gene polymorphism is associated with the lactation 
shift subsequent to delivery of male or female calf in 
domestic cattle. If allocation of resources according to sex 
of the calf allows optimizing the energy cost of lactation, it 
may result in positive selection of the G allele in domestic 
cattle breeding in extreme environmental conditions. In the 
same way, if the dramatic fall in milk production after birth 
of a heifer calf increases the probability of culling for the 
cows with the G/G genotype in animal husbandry, it may 
result in negative selection of the G allele during the 
subsequent lactations in dairy cattle herds. As regards the 
importance of our results to milk industry it seems 
reasonable to improve of cattle herds by selection in favor 
of A/A homozygotes. In addition to resistance to bovine 
leukemia, milk performance traits of cows with the A/A 
genotype are independent of from sex of newborn calf. 
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